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Abstract 



We explore the perturbative QCD corrections to the inclusive semilep- 
tonic decays of B mesons. In inclusive B decay processes the parton picture 
works well and the scaling feature manifests because the mass of b quark is 
larger than the QCD scale. Due to this property, the decay rate may be ex- 
pressed by a single structure function describing the light -cone distribution 
of b quark apart from the kinematic factor. We derive a q 2 evolution equation 
for the distribution function, which violates the scaling property through the 
Altarelli-Parisi type approach. The evolution equation is numerically solved 
and its phenomenological implications are also discussed. 
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I. INTRODUCTION 



Our understanding of inclusive decays of B mesons has made great progress over the last 
few years. Since the B meson is heavy, the momentum transfer of the decay is larger than the 
typical energy scale of hadronic bindings ~ Aqcd in most region of phase space and so the 
light-cone dynamics dominates. This enables us to formulate the inclusive decay precesses in 
an analogous manner to that of the deep inelastic scattering (DIS). The light-cone dominance 
allows us to express the commutator of two currents as bilocal operator of quark fields. The 
Fourier transform of the bilocal operator is related to the structure function of the decay 
which describe the momentum distribution of a b quark inside B mesons. Nonperturbative 
power corrections are encoded in the matrix elements of the bilocal operator and also in the 
distribution function. Because of asymptotic freedom of QCD, we expect that the B meson 
could be treated as a group of partons neglecting interactions between quarks following 
Feynman's parton picture U, which leads to the simple modellings of the inclusive decay 
processes 

With applying the operator product expansion (OPE) to the hadronic tensor, the leading 
twist contributions to the power corrections for the B decay processes have been extensively 
studied in the framework of the heavy quark effective theory (HQET) and the hadronic 
matrix elements are systematically parametrized in terms of 1/rrib expansions However 
the complete descriptions of the decay processes are not attained yet when the final state 
quark mass is small because of its singular behavior at the endpoint region where the OPE 
fails. Although the scaling property due to the heaviness of the b quark makes it possible to 
express the decay process by a single universal function, this function cannot be calculated 
from QCD at present. We still have unknown part of nonperturbative power corrections as 
well as we need to resum the leading twist contributions to the distribution function which 
becomes important at this region. Several authors have proposed the resummation of the 
leading twist contributions and smearing to obtain a well-defined shape function HTUHl^] . As 
a result we still lean on modelling the distribution function to take into account the decay 
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spectra at the endpoint region ||13|| . Therefore it is interesting to analyze phenomenological 



models in the framework of the HQET. The l/m^ structure of modelling the inclusive decays 
based on defining the mass of b quark as 

m b = (x)m B (1) 

are studied in the Refs. |||4| which leads to the parametrization of the nonperturbative 
contributions. Besides we can also find many attempts to the analysis of phenomenological 
models based on the HQET 

On the other hand there are perturbative corrections which have logarithmic singularities 
close to the endpoint at the quark level, which are regularized when integrating over the 
phase space and summing over ensembles of states. In phenomenological and theoretical 
ground, it is important to incorporate both perturbative and nonperturbative contributions 
at the endpoint region. In order to accomplish the task an evolution equation for the heavy 



quark distribution function has been studied in Refs. |I5 , I6[] . 

In the present paper, we attempt to derive an evolution equation for the distribution 
function in the intuitive way introduced by Altarelli and Parisi W71 to combine the per- 



turbative corrections into the parton model approach and solve the equation. This type of 
approaches avoids the formal tools of QCD but is completely equivalent. The phenomeno- 
logical implications of the evolution are investigated through the q 2 - and charged lepton 
energy spectra. 

This paper is organized as follows: In Sec. II, we briefly review the formalism for inclusive 



semileptonic decays of B meson. We show that the formulation by Jin and Paschos [|X9|j20 
is equivalent to the simple parton picture. The derivation of the evolution equation for the 
distribution function of the heavy quark is given by the method of Altarelli-Parisi in Sec. 
III. The equation is solved for the moments of the distribution function in Sec. IV and 
also solved numerically in Sec. V. We obtain the evolution of the distribution function and 
investigate the phenomenologies of the decay spectra. Finally we conclude our work in Sec. 
VI. 



II. INCLUSIVE SEMILEPTONIC B DECAYS 

We have the differential decay width for the process B — > X q lv 

J V _ G F\ V qb\ 2 j rrr^y ^Pl <^Pu (r> s 

d 2E B • a> (27r)32^(2vr) 3 2^' [) 

where is the leptonic tensor and the hadronic tensor W^ v is defined as 

W " V = \^l (2^^ (27r)4 ^ (Pa ~ g ~^ )(jB|Jt/f|X)(X|JljB) ' (3) 
which can be written in a more compact form 

= -\j d 4 x e iq - x (B\[J"{x), r\0)]\B) (4) 

by using a completeness relation \X)(X\ = 1 an d physical condition q > 0. The B meson 
state is covariantly normalized by 

(B(p')\B(p)) = 2E B (27r) 3 5 3 (p'-p). 

It would be helpful for the OPE that we relate W^ u to the absorptive part of the forward 
scattering amplitude 

= ~J d 4 x e iq - x (B\TJ^(x)r\0)\B) , (5) 

through the optical theorem W^ v = 2 ImT^. 

In most region of the phase space, we expect that the dominant contribution of the 
integration in Eq. (4) to comes from the region close to the light-cone, < x 2 < 1/q 2 



We write the commutator on the light-cone as [|T9| , pO 



(B\[J„(x), Jj(0)]\B) = 2{S mvP - ie^)(d a A q (x))(B\b{0)^(l - l5 )b(x)\B) , (6) 

where S mv p = g m g v p + g^g V a — g^uQa^ and A q (x) is the Pauli- Jordan function for the 
produced quark q defined as 

d 4 p 



iA q {x) = J 7^e-* p - x e(p )5(p 2 -ml 
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Note that Eq. (6) is Lorentz covariant and can be calculated in any reference frame. 

Light-cone dominance enables us to keep only the leading term in the light-cone expan- 
sion of the above reduced matrix element and we define the quark distribution function by 
the Fourier transformation of it: 

/(*) = / d{x-p B ) e-v B (B\b(0)p P Y(l - K)b(x)\B)\ x2=0 . (7) 

After integration over x, we obtain the hadronic tensor 

= 4n(S mi// 3 - ie mu p) J dz f(z) e(zp B0 - q ) 5((zp B - qf - m 2 q ) (zp B - q) a p P B . (8) 

We consider the hadronic tensor in the light of the parton model. Using the relation 
of decaying quark momentum pb = zp B and produced quark momentum p q = zp B — q, we 
rewrite the hadronic tensor in the form: 

W" v = I — f(z) W"" , (9) 

where 

is the spin sum of two quark currents like the leptonic tensor. We find that actually the 
hadronic tensor at the hadronic level is the convolution of the hadronic tensor at the quark 
level and it leads to the convolution equation of the differential decay rate 

dT(B -> X q lv) = J dz f(z) dT(b -> qlv) . (10) 

We find that Feynman's parton picture manifests in this approach. 
From Eq. (2), we obtain the triple differential decay rate 

dT G 2 F \V ub \ 2 q — E\ 



dE t dq 2 dq 47r 3 m B ^ q 2 + m 



f{x + ){2x + E l m B -q I )-{x + ^x„) , (11) 



where x± = [(q • p B ) ± y(q ■ p B ) 2 — m 2 B (q 2 — m 2 )]/m 2 B . The contribution of the second term 
dependent on x_ is negligible because the function f(x) should be sharply peaked. This is 



true after the evolution, which was checked numerically later. If we let m q be zero and take 
the B rest frame, finally we obtain the triple differential decay rate revealing the scaling 
behavior 

dT G 2 F \V ub \ 2 2 ^A(q ~E l )(2m b E l -q 2 ) 

-f(x+) , (12) 



dEidq 2 dqo 1927T 3 m B m b — go 

which effectively depends upon one scaling variable x + alone apart from kinematic factors. 
Here m b is as defined in Eq. (1). If we consider the perturbative QCD corrections, additional 
q 2 dependence comes in and the scaling behavior should be violated. We study the evolution 
of the distribution function of b quark with respect to q 2 in the next section. 

III. EVOLUTION EQUATION FOR THE DISTRIBUTION FUNCTION 

To incorporate the perturbative effects of QCD with our viewpoint of the parton model, 
it is very suggestive to apply the method of the Altarelli-Parisi equation fT7| , [T8|j since the 



inclusive semileptonic decays of heavy flavours are intimately related to DIS via channel 
crossing. This method will show close contact with the parton model approach to inclusive 
B decays at each step. We derive the evolution equation for the light-cone distribution 
function of b quark in this section. 

We decompose the B meson decay process into the subprocess as shown in Fig. 1 to 
write the amplitude []: 

M(B -> qXW) = M(B -> Xb)M(b -> qW)——l — . (13) 

We substitute the decomposed amplitude in Eq. (2), we have 



1 Note that we distinguish the produced quark q and other final states X in this section. Please see 
Fig.l (a). Namely we call the final states (momenta) the produced quark q (p q ), the rest hadronic 
degrees of freedom X (p x ) and lepton pair W (q). 
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dY(B -> qXW) = -?—\M(B -> X6)| 2 |A4(6 -> gWO| 2 



2^ B ' V 71 ' V 71 (#x + ^ - E B f (2E b y 

*(**)'*(*.-<-'* -^ \^n Ex{ 2%, ( ^ ■ (14) 

and comparing this equation with Eq. (10), we obtain 

f(z)dz = ^\M(B - Xo)| 2 1 -I ffi* , (15) 

where z = E b /E B . 

We consider the one-gluon emission from the b quark which brings out the leading cor- 
rection to the decay process of order of a s . Our strategy is to describe the effects of the 
one-gluon emission by the change of the distribution function. Let us consider b — > qlv 
process as a two-body decay with the final state quark q and the lepton pair (lu) = W. At 
the quark level, the differential decay rate is given by 

dYjb -> qW) 1 ,_2r 

= T 6{q - -m b ) = (16) 

ago z niB 

where = 2q Q /mB and q is taken to be a momentum of W boson (lepton pair). At the 
hadronic level, we have 

Thus we express the distribution function with an argument y + as the differential decay 
rate. The gluon emission yields the change of the distribution function 

Am = ^v (B ^ q xw + a) 

ZL o aqo 

As in the Eq. (13), we decompose the gluon emission process into three subprocesses 

M(B - qXW + g)=M(B^ Xb)M(b ^qW + — ^- , (19) 

h/X + &b — ZJ^ b 

M{b ^ qW + g)=M{b ^ h > g)M{h > _» qW)——± -JL , (20) 

as shown in Fig. 1 (b). The differential decay rates are given by 
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dT(B - qXW + g) = ^r\M{B - Xb)\ 2 \M(b - ^)| 2 |A^(6' - qlu)\ 2 

1111 

x 



(£ x + E b - E B f {2E b f (E g + E' b - E b y {2E' b f 
x (2vr) 4 5 4 (p B -px-q-Pq-Pg) 

d 3 p x d 3 Pq d 3 q d 3 Pg 
(2tt)32£ x {2nY2E q {2nf2q (2n) 3 2E g ' 1 ' 

and 

^ - w - - <w)\W(p« - , - p,) {2 % q(2 X» ' (22) 

Substituting Eq. (15) and (22) into Eq. (21), we obtain Af(x) from Eq. (18): 

A/M = W - x) ^)l 2 (^T^F (2lFpfe • <23> 

Now we need the amplitude for the b — > b'g subprocess to calculate A/(x) and it can 
be easily calculated. For simplicity, we take the limit \p b \ » m b to obtain the Lorentz 
invariant quantity \M.(b — > ^'g)] 2 . Let the components of the momenta of the b quarks 
before and after emitting a gluon and the emitted gluon be: 

rf = (yE B ,0,yE B ) , 

Pv = (v'Eb + kT, , (24) 

p^{(y-y')^ + ^^ B ,-^y-y')E B ) , 

where |k T | denotes the size of the transverse momentum of emitted gluon. We also let 
z = y'/y and obtain the amplitude 



v2 



= 87ra a C 2 (F) 1 + * |k T | 2 , (25) 
z(l — zy 

where e* fJl (p g ) is the polarization vector for gluon and C 2 (F) = Y, a T a T a is the casimir of 
SU(3) group. The factor 1/2 in the first line denotes spin average and 1/3 color average. 
The change of the momentum transfer arises from the gluon emission which sets the size of 
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|k T |. Hence we let <i|kr| 2 /|kr| 2 = dq 2 /q 2 and we have the contributions to the change of 
the quark distribution function: 

A/(:r, q 2 ) = f l dy dz f(yi tf) P{z) [6{yz _ x) _ 5(y _ x)] d(ln g2) ? (26) 

271 Jo 

where 

^) = C2{F) l -±^ (27) 
l — z 

represents the probability for a quark radiating a gluon such that the quark's momentum 
is reduced by a fraction z and has the same form as the non-singlet contribution of the 
splitting function of DIS. Note that we should have the negative contributions to Af(x) 
as well as the positive one. The negative contribution arises when the b quark originally 
has the momentum fraction x and the gluon emission reduces the momentum fraction into 
y' = zy. While the positive contribution term S(y' — x) implies production of a quark with 
momentum fraction x, the negative contribution term 5(y — x) implies elimination of a quark 
with momentum fraction x. 

Finally we write the evolution equation which is the same form as the non-singlet con- 
tribution of the Altarelli-Parisi equation: 

df(x,t) a a (t) f dy f x' 



dt 2,7c J y \y 

where t = \ In q 2 /q^ and 

:i + z 2 ) 3 



P(z) = C 2 (F) 



[l-z)+ 2 
with the notation of the function regularization 



+ -5(l-z) 



(29) 



o (l-z) + ~Jo d l-z ' m 
which denotes the regularization of the singularity due to soft gluon. This is introduced 
because this integral actually diverges as z — > I. We should also have included the vertex 
correction diagram for the infrared divergence cancellation when we squared the one-gluon 
emission diagram. 



IV. EVOLUTION OF MOMENTS 



We define the moments of the distribution function f(z) by 

M n (t)= l 1 dzz n f{z,t) , (31) 

JO 



which yields the evolution equation of moments 

dM n {t) a s (t) 



dt 2vr 

where 



M n (t)D n , (32) 



D n = f 1 dz z n P(z) . (33) 
Jo 

The evolution is caused by running of the strong coupling constant ct s (t) and Eq. (32) yields 



dM n D„ da x 



M n (11 - \nj) a s 
We obtain the solution for evolution of moments: 

M n (t) a s (t)\ o^FT) 



(34) 



(35) 



M„(0) \a s (0)J 

If we have an model distribution function at q 2 = q 2 as an initial condition, its moments 
M n (0) can be explicitly calculated and we obtain the evolved moments M n (t) directly. 

To show the result in terms of the HQET parameters, we consider the relations of mo- 
ments of the light-cone distribution function F(x) defined by the HQET and of our distri- 
bution function defined in Eq. (7). These functions are related by a change of variable 

-L/(, + ) = If (is) , (36) 

where 

ml + q 2 - 2m b q 

XB = ~ 2A(m t -,„) ' (37) 

If we define the moments of F(x B ) as 



10 



dx B x n B F(x B ) , (38) 

the first moment is independent of the scale according to the normalization of the function 

ao = 1 = M (t) , (39) 

which tells us the 6-number conservation corresponding to the moments sum rule of the DIS 
process. 

The second moment a\ is equal to be zero by the equation of motion of the HQET 
without running, which indicates the absence of the leading correction in l/m^ expansion. 
It holds when we take the pole mass of the b quark as the expansion parameter and is 
related to the definition of the b quark mass. In Ref. || it is shown that the parton model 
formulation can have this property if we define the b quark mass by the Eq. (1). Including 
running effects, however, Eq. (1) no more holds and the b quark mass should be redefined: 



m h = m B M x {t)[^%r\ ~ lUf , (40) 



which results in 



171 B li/r U\ I 1 ( 



""I"/ 



"M-X^l 1 -^. 1 1 (41) 

We can compute it using M 1 (t) and D x . 

The third moment a 2 is related to the average kinetic energy of the b quark inside B 
meson, fi 2 . So the running effects is expressed by the running of fj% and computed by 

f il(t) = 3m 2 (M 2 (t)-M 2 1 (t)) . (42) 



V. SOLVING THE EQUATION AND PHENOMENOLOGIES 

Since our evolution equation has the same form as the non-singlet part of the Altarelli- 
Parisi equation for the DIS, we use the package developed by Miyama, et al. pTI] for 



solving the Altarelli-Parisi equation numerically. This program adopted, so-called, the 
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brute-force method which is the simplest one in solving the integrodifferential equation. 
It divides the variable x and t into small steps and defines differentiation as df(x)/dx = 
(/(zm+i) - f(x m ))//\x m and integration as Jdxf(x) = E„ = i Al m/(im). Thus the evo- 
lution equation could be solved by repeated summations and subtractions. This method 
provides us the numerical accuracy as high as we want by spending running time enough if 
the initial distribution function f(x,ql) is given. At the low energy scale q 2 ~ ^qcd, the 
nonperturbative effects dominate and it is reasonable to take a modelling of the distribution 
function as an initial condition for the evolution at this scale. Here we use the fragmentation 
inspired form suggested by Lee and Kim || as the initial condition at q$ — 1 GeV 2 which 
improves Peterson's fragmentation function 



For numerical analysis, we took the values qlq = 0.085, tq = 0.004. In this case, the 
normalization constant is fixed as Nq = 0.0235 by the condition 



We solve the equation in the leading order of a s . 

The distribution functions for varying q 2 are shown in Fig.2. We can see that it becomes 
rapidly broader as q 2 increases, which suggest the QCD leading order resummation is large. 
We also show the plots of g 2 -dependences of the distribution function for different values 
of x in Fig. 3. We can understand these figures as follows: If x is large, the b quark has 
large momentum and would radiate more energetic gluons, which reduce the momentum of 
b quark. As a result, the low x region receives much positive contribution from the large x 
region and only a little negative contribution, while in the region around the peak, negative 
contribution dominates and positive contribution are very little. 

On the other hand, the observable such as E e distribution is not so much affected as 
shown in Fig. 4 because of the kinematical factor. We find that q 2 spectrum reflects the 
evolution effects of the distribution function more directly. It is recommended to measure 
the q 2 spectrum accurately for probing the evolution effects from perturbative QCD. 



f(x,q 2 = ql) 



N 



(43) 




2 




(44) 



12 



VI. CONCLUSION 



We derived and solved the evolution equation of the light-cone distribution function 
of a b quark inside B meson to combine the perturbative QCD corrections for inclusive 
semileptonic decays. In the viewpoint of parton model, the nonperturbative effects are 
encoded in the structure function of the decay process which discribes the distribution of b 
quark momentum. Since b quark is heavy enough, the momentum transfer of decay is larger 
than the QCD scale in most region of the phase space. It enables us to treat the b quark as 
a noninteracting parton. 

With the scaling property set in, the leading correction of perturbative QCD comes from 
the one gluon emission. The evolution of the distribution function is obtained by adding the 
one gluon emission diagram which violates the scaling behavior of inclusive decays. 

The distribution function obeys kinds of sum rules since the moments of the function can 
be directly related to the HQET parameters. If the mass of b quark is properly defined in 
the model framework, we can apply the l/m& expansion technique to the model study and 
that mass should be the pole mass. When the evolution effects comes into consideration, 
the pole mass is redefined to modify the equation of motion by the power corrections [JT6j . 



The fact that the second moment of the genuine QCD distribution function F(x) is equal 
to zero defines the b quark mass by Eq. (1) from the relation 

a 1 = ^( Ml -^)=0. (45) 
A \ m B J 

If ai is kept to be zero with running, the value of m^/i) = M 1 ( y u,)m fl changes and this 
rrib cannot be the pole mass. Therefore we should modify the definition of the pole mass 
to compensate the running effects as in Eq. (40) which indicates that a\ is no more zero 
with evolution effects. In our approach, we explicitly present that the definition of mass is 
changed by the running of ai(/i). 

We numerically solve the equation and the charged lepton energy and q 2 spectra are 
also shown. The shapes of decay spectra are changed by incorporating the perturbative 
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corrections, while charged lepton energy spectrum is less sensitive. The decay rate decreases 
with the corrections and we conclude that the perturbative corrections should be considered 
in practically extracting the CKM parameters. 
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FIGURES 

FIG. 1. (a) Diagram of inclusive semileptonic decay of B meson, 
(b) One gluon emission diagram in Inclusive semileptonic decay of B meson. 

FIG. 2. The light-cone distribution functions of b quark for q 2 = 1.0 (solid), 5.0 (dashed), 25.0 
(dash-dotted) GeV 2 . The curve for q 2 = 1.0 GeV 2 is the initial distribution. 

FIG. 3. The evolutions of the light-cone distribution functions for fixed x. From up to down 
x = 0.9, 0.8, 0.7, 0.5, 0.1. 

FIG. 4. The g 2 -spectra of inclusive semileptonic B decays: The solid (dashed) line denotes 
the case without (with) the QCD evolution effect. Only the q 2 > 1.0 GeV 2 region from which the 
distribution function is evolved is shown in the figure. 

FIG. 5. The electron energy spectra of inclusive semileptonic B decays: The solid (dashed) 
line denotes the case without (with) the QCD evolution effect. 
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